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Abstract 

We give a detailed and rigorous picture of the mesoscopic con- 
ductance fluctuations in the deep insulating regime (DIR) within the 
Nguyen, Spivak and Shklovskii model including spin-orbit coupling 
(SO). Without SO, we find that fluctuations of the log-conductance 
are persistent above a saturation field Bg, where one has that the log- 
conductance is approximately a stationary random process. In con- 
trast, in the SO case the saturation field Bg is negligible and the sta- 
tionarity is well realized. We find non-vanishing disorder fluctuations 
of the field average of the log-conductance as a quantitative measure 
of the lack of ergodicity in the mean square sense. To this fact, a weak 
decaying behavior of the correlation function, even weaker in the case 
of SO is established on the relevant field scale of the model, one flux 
quantum per plaquette. This finding corroborate the behavior of the 
fluctuations of the field average of the log-conductance and permit us 
to invoke Slutski's theorem to conclude that the whole stochastic pro- 
cess defined by the log-conductance is non-ergodic in the mean square 
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sense in both cases. As a consequence the commonly used criterion to 
test the ergodicity based on the equivalence of the variance in disorder 
and the variance in the field is not fulfilled. Using the replica approach, 
we derive the weak localization analogs of the 'cooperon and diffuson 
which permits us to analyze in qualitative form the decaying behavior 
of the correlation function. Our predictions agree qualitatively and 
semi-quantitative with experiments in the DIR. 
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1 Introduction 



The nature of Fluctuations in both the metalhc state and in disordered 

insulators 0, [§, ||, has been a matter of interest for both theoretical and 
experimental studies. Whereas in the metallic regime the basic aspects of 
fluctuations have been elucidated, in the regime of hopping transport the 
nature of fluctuations is still an open field. The deep insulating regime, DIR, 
where transport occurs via variable range hopping (VRH), is defined as the 
regime where the localization length is the smallest scale compared to the 
elastic mean free path and hopping lengths, i.e., C, < i < t respectively [^]. 
Coherence effects are possible in this regime because phase breaking events 
occur at the hopping length[^, which is larger than £. Important signatures 
of quantum interference in disordered insulators are the classic magneto- 
fingerprints, or reproducible fluctuations in the conductance with magnetic 
field, and a low field positive magneto-conductance. 

An important property of mesoscopic conductance fluctuations in the 
metallic phase is their ergodicity. At the mesoscopic level, the sample size 
is less than thermal diffusion length or the dephasing length, whatever is 
shorter, such that sample to sample fluctuations are visible and the system 
does not self-average. Although it was not rigorously proven, the ergodic 
hypothesis was meant as the ability of the magnetic field ( or energy) to 
induce conductance fluctuations equivalent to sample to sample fluctuations 
(Lee-Stone criterion) [Q. In contrast, experimental results show that log- 
conductance mesoscopic fluctuations in the DIR without spin-orbit scattering 
are not ergodic in the Lee-Stone sense [^, ||, ||, i.e., the variance over samples 
is larger than the variance over field. Such samples involve hopping lengths 
that are, at most, 6 to 10 times the localization length. Precise measurements 
of Ladieu et al[Q and Orlov et al0 have shown that a) field fluctuations do 
not decorrelate disorder fluctuations, b) field fluctuations do not change the 
identity of the hop, c) the field average of the variance over the samples is 
larger than the sample average of the variance over the field and d) there 
exist a decorrelation field Be defined by the field correlation function, which 
defines a equivalent new sample. 

The question of the ergodic nature of fluctuations with and without SO 
has not been addressed, to our knowledge, from the theoretical side. Our plan 
of this work is first to address the problem of fluctuations and the question 
of ergodicity in DIR within the NSS model. We undertake this task through 
the verification of concepts concerning ergodicity, which we first define with 
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mathematical rigor in the next mostly technical section. Then in section 
^ we explain the NSS model and define the random processes we want to 
analyze. In section^ we carry out the program described in section ^ and 
verify the non-ergodic behavior of fiuctuations. In section § we define the 
main theoretical objects of this work, the cooperon and diffuson analogs 
of weak localization theory, with the help of which, we can explain to some 
degree the decaying behavior of the correlation function. Finally we conclude 
by discussing and comparing our results with experiments. 



2 Ergodicity of Transport Fluctuations 

In the preceding section one we introduced the question of ergodicity of 
fiuctuations. Here we introduce the mathematical concepts that will permit 
us to establish the ergodic nature of the fiuctuations. 

Given a physical quantity F{T-C, B) depending on the disordered Hamil- 
tonian Ti and magnetic field S, we denote by F{7i, B) the sample to sample 
average, or disorder average, and by {F{n,B)) = AB-^ JsUB F{n,B) 
the field average for a given sample or disorder realization. In order to 
be able to estimate the sample average from the field average of a given 
sample the following conditions must be satisfied: a) lim^^^oo o"mss(-S/) = 

[F{n,B) - {F{n,B))f and b) F{n,B) = {F{n,B)). The verification 
of both conditions is known as ergodicity in the mean square sense (mss), 
or the random function -F(7i, B) is said to be ergodic in the mean-square 
limitp. The condition F(Ti.,B) = {F(7i,B)) is a measure of global sta- 
tionarity. which means that these averages are independent of B. One can 
cast conditions a) and b) into a single statement on the disorder fiuctuations 
of the field average,i.e., limB^_,oo cTmss (5/) = lims^^oo Varrf((F(7^, 5))) = 



{{F{T-C, B)) — {F{T-C, B))Y 0, where Var^^ means variance over disorderjrU 



This property implies that for one realization of disorder there are enough 
equivalent samples within the magnetic scale, such that the average in the 
field, with regard to disorder, does not depend in statistical sense on the 
particular realization. This means that (F(7i,i?))^ N{F{'H, B)), where 
is the number of realizations, and a sharp distribution of {F{T-l,B)) over 
disorder holds. 

One could also ask for the possibility of making estimates of Vaid{F{T-C, B)) 
from Va.TB{F{l-C, B)), ( here Vars(F(7^, i?)) means the variance over the 
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field B of F{T-C,B)), or more generally, try to make an estimate of an- 
other function of the basic process F(T-C,B), g{F(H, B)). The necessary 
and sufficient conditions such that one can estimate g{F{T-C, B)) from one 
realization of disorder with the field average {g{F{TC, B))){ the so called 
law of strong numbers), is given by Slutski's theorem. One can write: 
amss{Bf)=\imB,^oo ^ Ib! dB{Bf - B)C{ g{F{n, B))), with: 

C{g{F{n, B, AB))) = A[g{F{n, B + AB)) x A[g{F{n, B)) (1) 

where C{g{F{n, B, AB))) is the correlation fun ction for g{F{n,B ,AB)) 
and A[g{Fln, B + AB)) = [giF{n, B + AB)) - g{F{n,B + AB)) ]. One 
can easily realize that a strong decaying behavior of the correlation function 
with correlation lenght t <^ Bf will be a sufficient condition for ergodicity 
in the mss, i.e. amss{Bf) =^ 0. In fact, it is also a necessary condition. 

Usually, for the application of the theorem stationarity in the wide sense 
of F(H,B) can be assumed ,i.e., {F(TC,B)) does not depend on B and 
C{g{F(H, B, AB))) should depends only on A_B|0. In this work we are 
interested in testing the ergodicity in two cases: g{X) = X for the average 
and g{X) = X^ - X for the variance. The first case corresponds to the 
usual meaning of ergodicity in statistical mechanics. The second case, usu- 
ally named the Lee-Stone criterion, refers to the equivalence of sample to 
sample fiuctuations and magnetic field fiuctuations. [p]] P||1T2[| . 



3 The NSS Model 

We examine now the fiuctuations in the DIR in two and three dimensions 
and define the random process -F(7i, B). This is obtained from the Nguyen, 
Spivak and Shklovskii]^, ^ model (NSS). The NSS model's crucial insight is 
that coherence is maintained within a Mott hopping length, where the con- 
ductance is a sum of coherent forward directed Feynman paths which interfere 
which each other. The NSS model describes the quantum behavior of the 
critical (bottleneck) hop in the Miller- Abrahams network |]13||. The existence 



of many randomly oriented critical hops tend to average the macroscopic 
conductance, eliminating fiuctuations|]TB[. Here, we focus on the low temper- 



ature regime where critical hops do not trivially self average[3]||17||, i.e., the 
percolation correlation length is such that = ^{T^IT)^"'^^^/^^^^'' ~ L, 
where v is the percolation correlation length exponent, D is the spatial di- 
mension and To a disorder parameter. This is the mesoscopic regime [Q. We 
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will first find the fields above which the process F(T-C, B) can be considered 
stationary and test the ergodicity criterion in the mean square limit. We 
then calculate the correlations functions and their decaying behavior. Invok- 
ing Slutski's theorem, we confirm non-ergodic behavior as established by the 
finiteness of the mean square criterion, (lim^^oo Vard((F('H, 5))) 7^ 0). We 
then proceed to test the ergodicity of fiuctuations. Non-ergodicity of fiuctua- 
tions is implied by the weakly decaying behavior of the correlation function. 
We find this consistently. Here, we find good quantitative agreement with 
the measurements of Orlov et al p. We derive the cooperon and diffuson 
analogs of the weak localization theory for DIR with the help of which the 
non ergodic behavior can be explained. Furthermore , the predictions for the 
case with SO are made. 

In the two dimensional NSS model, impurities are placed on the sites of 
a lattice of main diagonal length t (the hopping length, t = ^{To/TY^^^^^\ 
Mott's law). We apply a magnetic field B, perpendicular to the plane, chang- 
ing only the phases of the electron paths. The overall tunneling amplitude is 
computed by summing all forward directed paths between two diagonally op- 
posed points, each contributing an appropriate quantum mechanical complex 
2x2 matrix weight given by the Hamiltonian: 

i <ij>a,crr 

where is the site energy, and Vij^a,ai represents the nearest neighbor cou- 
plings or transfer terms which includes a randomly chosen SU(2) matrix 
describing a spin rotation due to strong SO scattering. Within the NSS 
model, we choose site energies to be = ±W with equal probability |^ |]24 



Without SO the coupling terms are diagonal in spin space Vij = V, and the 
Green's function between the initial and final site is given by 

/ \/ \ t directed 

(^|G(E)|/)= (-) J(i?,t); J{B,t)= E [Uv^r,eM, (3) 

r i-p/ 

where is the phase gained through path ip' due to the magnetic vector 
potential, T' represents all directed paths that go from i to f through the 



lattice and T]i = sign (ej) = ±ip3|. In the presence of spinorbit scattering 



the Green's function is the 2x2 matrix: 

directed 
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The Green's function consist of sum of terms, one for each path, each being a 
product of random numbers and random SU(2) matrices, and a determinis- 
tic disorder independent phase factor from the magnetic vector potential p3 



The complex function J'^{B,t){ a complex matrix function in the presence 
of SO) contains the interference information including correlations due to 
crossing of paths, and the factor {V/Wy is the leading contribution to 
the exponential decay of the localized wavef unction. We use the trans- 
fer matrix approach in order to compute J{B,t), exactly, for each real- 
ization of disorder ||2^. Our random processes in question represent the 
log- conductance. They are F(T-C,B) = \n{J^{B,t)J{B,t)) and with SO 
F{n,B) = ln(/(S,t)) , where I{B,t) = l/2Tr{.P{B,t)J{B,t)) H. In his 
work we measure the magnetic field B or changes in magnetic field AB in 
flux units (po/^"^, where 0o is the flux quantum. 



4 Fluctuations and Ergodicity 

Figure |l| shows typical fluctuations of the log- conductance as a function of 
the sample and the magnetic field. Without SO , the figure clearly shows 
that the average of ln{J{B,t)) dominates the fluctuations, i.e., the average 
behavior is visible for a single sample. As the average ln{J{B,t)) first in- 
creases proportional to B, crossing over to a slower growth as B^^'^ dictated 
by the magnetic length AB < B^ = 7ich/{^^/^e t^/'^) [g§ ^, |?7|, the process 
is not stationary. However, in the latter regime of slow growth, one finds 
a field above which the process can be considered as essentially stationary 
in the same fashion as in the metallic regime. In this regime, while the 
log- conductance tends to saturate, the fluctuations persist as in mesoscopic 
fluctuation theory in metals[|^]. Furthermore, we note that the average be- 
havior is periodic in half the flux quantum 0o per f' (only one half of a 
period is shown). This periodicity reveals an average field coupling to 25 1|^ 
which has been demonstrated theoretically||31||. In three dimensions, the 
fluctuations are appreciably larger than the average behavior. Once more, 
persistent fluctuations beyond the average conductance saturation field are 
observed. The existence of such persistent fluctuations were first surmised 
by Sivan et al[|16], ^ and Zhao et al[0. With SO , there is no tendency to 



build an average, and there are in general soft changes in the fluctuations in 
marked contrast with the sharp changes without SO. This peculiarity will be 
explained later. One marked feature of figure l(with and without SO), are 
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the fact that disorder fluctuations do not decorrelate the field fluctuations, 
which is the reason of the remarkable similarity with fig. 2 of the work of 0. 
The key point in both figures is that the the fluctuations do no decorrelate 
at the scale of fields shown in the figures, suggesting non-ergodic behavior. 

Following the concepts of section two, we analyze the ergodicity of the 
log-conductance fluctuations more carefully. To achieve this we have to 
check two points: first if stationarity is reasonably fulfilled and second, we 
have to verify the ergodicity condition and the decaying behavior of the 
correlation function. Fig.||a shows the quantities Varrf((F(7i, (stars) 

and [F{n,B) - {F{n,B))]^{ diamonds) in the absence of SO {F{n, B) = 
ln\J{B,t)\'^). The field averaging interval is AB = [Bi,Bf]. In this figure 
as the field Bi is increased from zero to 0.09 in [0o/^^] units, the diamonds 
move downward overlapping the stars sooner for smaller t. Therefore, the 
field above which the process can be considered quasi stationary increases 
with the hopping length t such that both quantities tend to coincide for 
larger Bi . This feature is of importance for calculating other quantities, as 
any question on ergodicity presumes at least quasi-stationarity. In the case 
of SO iF{n,B) = ln|J(5,t)n, Fig.||b shows effective stationarity of the 
process \n{I{B,t)), essentially independent of Bt and t. Yaxd{{F{7i, B))) 
does not tend to zero with increasing Bj; on the contrary, it increases with 
a power of t|31], which is an indicator of non-ergodic behavior as there is no 
self- averaging as one increases the hopping length t( see Kramer and Mack- 
innon 0). The last result establishes non-ergodic behavior in the mean 
square sense with and without SO. To further substantiate this result we 
calculate the correlation functions. Fig. |^ shows the correlation function 
for three values of the hopping length t and some values for Bi. One sees 
a very weakly decaying behavior on the physical field scale( 0o/^^) and a 
tendency to decay faster for bigger Bi , indicating that when the process 
becomes quasi stationary there is a tendency to a faster decaying correlation 
function. Fig. ^3 with SO shows that the correlations function depends 
essentially only on AB, and the decaying is even weaker with a functional 
form depending on t. Now, the basic argument against ergodicity is that the 
decaying behavior of the correlation functions is such that it is not possible 
to construct enough ensembles from the field fluctuations data and there- 
fore non-ergodic behavior is established. To illustrate this point, we write 
(^mss{Bf) = ^Jsf dB[{Bf - B)C{ F{n,B))]. In order to have enough 
emsembles in the field scale within the validity of the model {Bf < 0o/^^ = 1 
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in our units) , there should be a decorrelation field << 1, such that large 
number of samples could be defined, so that one would get <Jmss{Bf) ^ 0- 
could be satisfied. This, on the other hand, implies that the condition for 
ergodicity of the variance, which in general requires a stronger decaying be- 
havior of the correlation function for process F{T-C, B) than the condition 
required for ergodicity in the mean square sense, is not at all fulfilled 0. 
Therefore, one should expect that the Lee and Stone criterion of ergodicity 
on the relative magnitude of the field and sample fiuctuations is not realized. 
To check this point one has to compare the magnitude of the variance in field 
and sample to sample fiuctuations. The idea of further averaging over dis- 
order and over the field, respectively, is the same as in statistical mechanics 
using different initial conditions to improve the statistics. Figure ^ shows the 
averages 

(Var,(ln|J(i?,t|)) = ((In 1 7(5, t| - In | J(5, (5) 



VarB(ln|J(i?,t|) = ((ln|J(5,t)| - (In | J(5, (6) 

Without SO we find the characteristic Bi dependence shown in fig.|^ shows 
up as a crossing of both types of averages. There is again a clear tendency 
to saturate as Bi increases whereas the difference of the saturated averages 
widens with increasing t. This tendency is clearly seen in the case with SO 
where there is a much weaker dependence on Bi as expected from fig.^] and 
the crossing observed in the previous case is absent. 

5 The Correlation function: the Cooperon 
and DifFuson 

In this section we develop the concepts of the cooperon and diffuson in the 
context of strong localization. These objects are used to explain qualitatively 
and semi-quantitatively non-ergodic behavior in the mean square sense and 
the relative magnitude of the field and sample fiuctuations found in experi- 
ments. For this purpose, the following relation is straightforwardly derived 

Var4F(5+A5, t)+F{B, t)] = VaTd[F{B+AB, t)]+VaTd[F{B, t)]+2C{F{n, B, AB, t)) 

(7) 
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This relation is valid for both processes In J{B,t) and lnJ(i?)|]T^. For 
the sake of clarity, the composite process inside the brackets on the left 
side of eqn.0 is denoted by the process P{B,AB,t) = [ln\J{B + Ai?)p + 
ln|J(E)|2] and with SO by Pspinor{B, AB,t) = [ln/(5 + AB) + \nI{B)], 
such that the left hand side of eqn. ^ is in each case Var d[P {B, AB,t)] and 
VaidlPspinoriB , AB ,t)], respectively. These two functions are shown in fig. | 
as a function of AB in (po/i'^ units and three values of t (30, 100, 300) from 
the bottom to the top, respectively). For Bi shown in the figure and the 
scales of AB used, three things deserve explanations: first for large enough 
AB, the above mentioned functions show almost no dependence on Bi and 
AB at fixed t; secondly, the ratio Varrf[P(_B, AB, t)]/ VaXdlPspinoriB, AB, t)] 
is around two, for big enough t and AB which is a landmark of the symme- 
try changing from unitary to simpletic; Thirdly, one observes a very rapid 
decaying on the dependence on B, which can be traced to the saturation 
behavior of the cooperon as we will see. 

Recall that on the metallic side, an analogous behavior has been de- 
scribed which identifies two fundamental contributions to the field effect: the 
cooperon and the diffuson[|l|0. These contributions can be distinguished by 
the way they enclose the magnetic flux; while the cooperon is sensitive to 
{2B + AB) , the diffuson only responds to field changes AB. In the insulating 
regime, a mechanism similar to the cooperon which saturates is associated 
with a positive magneto-conductance (MC). This has been observed as a 
general effect |T5|, ^ |T^. A semi quantitative explanation for the behav- 
ior of the functions Var^ [P(i?, Ai?, t)] and VsndlPspinoriB, AB,t)] and their 
ratio can be found with the help of the cooperon and diffuson analogs. To 
achieve this goal, we consider the moments of process of P{B, AB,t) and 
Pspinor{B, AB,t){ further below it will become clear why). In the former 
case they are given by [J*{B + AB)J{B + AB)J*{B)J{B)Y. Recall from 
equations |^ and |^ that this product can be visualized as a set of n paths, 
each one defined by the respective term in the product. These paths even- 
tually intersect each other at some values of their length. In order to have 
nonzero contributions after disorder average, the paths must pair up, as a 
consequence of th chosen distribution of the energies H. Neutral paths 



(field independent) are formed by pairing J* and J at the same field (phase 
cancels). On the other hand, charged paths (field sensible) are formed by 
pairing either J*{B + AB) and J{B) or J*{B + AB) and J*{B). In the ab- 
sence of paired path intersections, self interference kills charged paths (their 
contribution decays exponentially fast). Nevertheless, if intersections are con- 
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sidered, one can have path exchanges for short distances, yielding a magnetic 
field coupling which is analogous to the magneto-conductance, the source of 
the initial decaying of the correlation function. There are three possible 
diagrams at a paired path crossing, two of which are depicted in fig. 7. With- 
out SO the spin indexes can be ignored, one obtains [^]: a) one partner 
from J*{B + AB) pairs with one from J*{B) while one from J{B + AB) 
and one from J{B) follow a different path. Such a combination encloses 
{2B + AB) and is therefore called cooperon-like. b) One partner is taken 
from J*{B + AB) and the other from J{B) on the same path, while one from 
J{B + AB) and J*{B) follow another. Such a combination encloses only AB 
and is called dijjus on-like. Finally, one can have combination c) where one 
partner comes from J*{B + AB) and the other from J{B + AB), leaving 
J*{B) and J{B) to pair up. The latter combination is called uncharged and 
encloses no field. Note that all previous cases satisfy overall neutrality so 
that the contributions are real as expected. The contribution of the replica 
cooperon and diffuson are the same at zero field and there is an additional 
contribution from the uncharged diagram. Further progress is achieved using 
the replica argument. 



The replica-moment argument p9| , maps the n-th moment problem onto 



the problem of 2n bosons with contact interaction. These interactions renor- 
malize due to the diagrams above, making path interactions field dependent. 
The 2n boson system can be solved using the Bethe ansatz and has ground 
state energy eo = ln4"' + p{B, AB)n{n^ — 1), where p{B, AB) is a function of 
5, AB, such that one has [J*{B + AB)J{B + AB)J*{B)J{B)Y =A{n, B, AB) exp(/n4"+ 
p{B, AB)n{n^ — l)t ) valid at fixed n asymptotically for t —>■ oo. On the other 
hand, the n-th moment can be expressed as a cumulant expansion valid 
at fixed t asymptotically forra -> 0, [J*{B + AB)J{B + AB)J*{B)J{B)Y 
= exp{E fCi[P{B, AB, t]}, where Ci[P{B, AB, t] are the cumulants of pro- 
cess P . The subtleties concerning both limits have been discussed by Kardar 
||29|| , who finds a nonextensive correction subleading term proportional to 
^2/3_ therefore obtain : 



Var,[P(5, AB, t)] = {pcoop{2B + AB) + pa.ff{AB))t^/^ + In A{B, AB) (8) 

Here p{B, AB) = pcoopi'^B + AB) + pdiff^AB), we have separated the path 
interaction in terms of the cooperon and diffuson contributions. We check 
numerically this important prediction and fig. 6 shows the expected scaling 
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with an exponent near |. Beyond the saturation field of the average log- 
conductance, the cooperon term on the right hand side saturates (the same 
for the variance on the right of eqn. |^ which depend on B) and the correlation 
function only depends on AB. This behavior is summarized in Fi^. 

In the case of SO, it can be shown |]2^, that the only non-zero paired 
averages are UapU'^p = |, t^n^u ^ h ^U^h = ' ^^^^ averaging 
brings a factor of (|)^ and forces the neutral paths to have parallel spins 
while the spin of the two partners of charged paths must be antiparallel. As 
a consequence, one finds the cooperon diagrams cancel in pairs as concluded 
for the case of the magneto- conductance [0] (fig. 7) such that no exponential 
corrections to the conductance occurs due to the cooperon. For the diffuson 
there are only two combinations possible for incoming spin indexes(all up and 
all down), such that in this case VaTd[PspinoriB, AB,t)] ={p^£jy^{AB))t^^^ 
-|-lny4(i?, AB). This nice result explains why there are fiuctuations with SO, 
even if there are not exponential corrections to the magneto-conductance. We 
find again numerically the predicted t'^^^ scaling (see fig.6). For large AB and 
t one finds the ratio of the variances approaches ~ pdiff{AB))/p^^jy^{AB) = 
2 (i.e., 1/(2 times (1/2)^)), in agreement with the numerical results. Now 
from eqn. |^ one obtains : 



C{B,AB,t) = ^{[(p,oop(25+AS)+prf,;^(AS))-(p^e(S)+p„,(S+A5))]t2/3+^(5,A5)} 

(9) 

where 5(5, A5) = In A{B, AB) -In A (B+AB) -In A (AB), are logarithmic 
corrections from the prefactors. pmc{B) defines the magneto-conductance, 
[J*iB)J{B)]- =exp{E fCSnJ*iB)J{B)]} =A{n, B, AB) exp{A\nn+pmc{B)n{n^- 
l)t) 0. 

From eqn.^ apart from the predicted t"^^'^ scaling, one can gain qualitative 
and quantitative understanding of the decaying behavior of the correlation 
function for small AB. With the explicit field dependence of the cooperon 
and diffuson given by: 

Pfiuctuauons = (2 Vs) (cos( (25 + AB)/B'^) + cos(A5/5^) + l)p{B = 0)t2/3 

(10) 

and 

Pmc = (l/3)[(2+cos(5/S,)) + (l/3)(2+cos((5+A5)/5,)]p(B = 0)^^/3 (11) 
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one obtains: 



C{B, AB, t) = l[8 cos((5 + AS) / B[) cos{B / B[) - 2 cos((25 + AB)/2B^) cos(A5/2fie)]Var4F(0, 
6 

(12) 

Here, we have ignored the logarithmic corrections from eqn.^. For i? = 
and = one has C{B = 0,AB = 0, t) = Var4F(0, t], with Var4F(0, t] = 
p{B = 0)t^/^, in accordance with eqn.|^ . Now, for small AB and i? ~ 0, one 
can assume that B^ = i?c |2^: 



CiB ^ 0, AB, t) = [l- ^(Afi/5e)']Var4F(0, t]t^/' (13) 

This behavior is is qualitatively seen even for i? 7^ in fig.^, where for given 
B one observes a faster decay with increasing t{ recall that B^ = vr/t^/^^^j^^ 
This is even more clearly to see with SO. In that case one has: 

a,_(i?, AB, t) = ^[(pI}7(Ai?) - 2pr"°1t'/' (14) 



Csp^nor{B,AB,t) = i([(2V2)(cos(Ai?/i?:) + 1)] - 2)pr"°'^t2/3 ^15) 
again an expansion in small AB : 

C^j^noriB, AB, t) = {l- ]^{AB/B,)^)ptre'^ (16) 

On the other hand, without SO, numerical calculations (fig.^, shows 
a very slow decaying, such that for the field range less than (po/^'^i there 
are no enough ensembles to form, such that even if the field correlation 
function decays, and a decorrelation field can be defined, it is not possible to 
define a large enough number of samples such that averaging over the field, 
could be equivalent to averaging over disorder (samples), a matter that has 
caused some confusion. With SO the decaying of the correlation function is 
dominated by an even much bigger field scale unknown to us. 

Although the replica results are valid quantitatively for 1 + 1 dimensions, 
nevertheless, 2 + 1 dimensions also has a bound state, although weaker, with 
a smaller positive MC. Therefore, our results apply qualitatively to three 
dimensional hops[^ . 
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6 Discussion and Conclusions 



Summarizing the results of this work, we have found that the log-conductance 
in the DIR of mesoscopic samples is non-ergodic in the mean square sense 
and the fluctuations of the log- conductance are non-ergodic in the sense 
that sample to sample fluctuations are larger than magnetic field fluctua- 
tions. Without SO due to the ever increasing magneto-conductance with the 
field ||2^, only quasi stationarity can be achieved for fields larger than 
a certain hopp dependent field Bs{t), whereas with SO due to the small 



magneto-conductance and its rapid saturation ||2^, Bs{t) is negligible, so 
that the process in that case can be regarded as stationary independent of 
t. The field Bs{t) is obtained numerically by varying Bi until one observes 
that amss{Bf) ~ Vard((F(7^, -B))) . This different behavior with and with- 
out SO turns out to have a remarkable influence on the behavior of the two 
variances that define the sample to sample or disorder fluctuations and the 
field fluctuations, as their values depend on B^. It is only when the averaging 
interval is taken as AB = [Bs{t), 5/], i.e., when quasi-stationarity is achieved, 
that the proper comparison of the variance can be made. In experiments the 
hopp values were small, and therefore Bs{t), was small. However, in order to 
compare the magnitude of the fluctuations defined by eqn.|^ and eqn.^ with 
experiments one must be aware that their values are sensitive to [Bs[t),Bf], 
an important point indicated by our results. The experiments of were 
done without SO; therefore, in order to make a a rigorous comparison with 
our results we have to know the experimentally taken interval [Bs{t),Bf] for 
the evaluation of this quantities. For the small values of t, according to our 
results, the influence of B^ is small and one should expect Bs{t) ~ and 
indeed we find the ratio of eqn.^ and eqn.^ predicts a value similar to their 
experiments. On the other hand, the scaling of the disorder fluctuations with 
the hopp length t, as shown in fig. 4, agrees with the experimental values ob- 
tained by Orlov et tal[Q(see fig. 45 in the review of Kramer and Mackinnon). 
It is interesting to observe that both the variance with disorder and the one 
with field appears to have the same functional form. A more careful com- 
parison can be done by taking into account all predictions of this work. In 
the case with SO, experimental studies of fluctuations are lacking. To our 
knowledge the only work where field fluctuations in the strongly insulating 
regime have been observed in samples with SO is the work of Hernandez and 
Sanquer p2|. Here our predictions give a clear qualitative way of differen- 
tiating the cases with and without SO and highlight a way of evaluation of 
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experimental data. 

From the theoretical point of view, we have derived two important ob- 
jects: the cooperon and the diffuson which are the weak localization analogs 
in the DIR. They allow us to explain some features of the correlation func- 
tion like the qualitative behavior for small Ai?, and numerically the weak 
decaying behavior on the scale (0o/^^), that according to the Slutski's theo- 
rem, is responsible for the non-ergodic behavior as we have found. In order 
to experimentally observe persistent diffuson fluctuations, one has to explore 
a range of parameters so that there is a saturation in the average behavior 
while the wave function shrinkage ||13|| is still a negligible effect. This range 
can be defined by the condition Be < h/ [easN^^^] = B^rb, where is the 
Bohr radius. Borb is the scale for the orbital shrinkage to be important [Q, 
i.e., when the cyclotron radius becomes of the order of the mean free path i. 
These conditions have been met in ref. |0 and ||^. Furthermore, according 
to references[||, ^ || the magnetic field cannot induce geometric fluctuations 
due to changes in the identity of the hop ||3CI|| . This finding holds in both two 
and three dimensions, and therefore we expect the insulating cooperon and 
diffuson fluctuations should be seen experimentally also in three dimensions. 
In experiments , one should be aware about the condition of mesoscopic 
sample, discussed in section ^ Otherwise, trivial self-averaging of spatially 
different oriented hops can wash out the possibility to observe fluctuations. 

There will be a non-ergodic to ergodic behavior in the case when one 
relaxes the condition of DIR with i << < t such that there are many 
impurities within C,- In this case there will be a diffusing behavior within 
the length scale C, so that two overlapping random processes are at work, 
one ergodic for the diffusing scale ^ and the other non-ergodic for the larger 
scale t. This question is under study |jl2|. Finally, we want to stress that 
the fundamental point of considering correlations in the random processes 
F(H,t) due to the crossing of paths can not be relaxed. This approach 
permitted us to use the full power of the replica theory, which again helped to 
semi-quantitative and qualitative predictions. Theories like the independent 
path approximation, where such crossing is neglected, miss the very crucial 
objects defined by the cooperon and the diffusson. 
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Figure 1: Field Conductance fluctuations of ln\J{B,t)\ and ln\I{B,t)\ for 
different realizations of disorder or sample number and t = 15 
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Figure 2: Quasi Stationary behavior without SO. Determining Bs{t) for 
t=30,50, 100, 150,300. Diamonds converge to stars symbols as the field Bi 
grows from zero to 0.09 in 0.01 steps. In the case with SO diamonds and 
stars fall together independent of and t indicating stationarity with Bg ^ 

Figure 3: Correlations functions C{B, AB, t) and Cspinor{B. AB. t) as a func- 
tion of AB{(j)/(j)o) for a given value of B from B = 0.0 to B = 0.09. For 
each value of B, there are three curves corresponding to three values of 
i(30, 100, 300). One sees a tendency to faster decaying with increasing t . 
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Figure 4: The figure depicts the variance as defined by Eq.| and Eq.|^ as a 
function of t. The effect of changing the values of Bi from to 0.08 is clearly 
seen until saturation is attained. Notice that there is a crossing due to a 
stronger with t increasing dependence on Bi{t) on eqn.^ than in eqn.^. In 
the case with SO, both quantities have a weak dependence on Bi and the 
crossing efect is absent. 
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Figure 5: The variance in disorder of P{B,dB,t) and PspinoriB,dB,t) as a 
function of dB{<p/<po) for t = 30, 100, 300 and B = Bi = 0.0- 0.09. 
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Figure 6: The figure shows the expected t^/^ dependence from eqn.§. Each 
hne corresponds to a different AB. The upper hne corresponds to AB = 0. 
With increasing AB the hues move downward and saturate for big enough 
AB. One observes small modulations which are interpreted as deriving from 
the normalization factors in the Bethe Anzatz. 
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Figure 7: The Cooperon and Diffuson without SO 
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Figure 8: The Diffuson and the Cooperon with SO 
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